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ABST RACT 

We use iDiaferid (|l999T >'s caustic method to identify the member stars of five dwarf 
spheroidal (dSph) galaxies of the Milky Way, the smallest dark matter (DM) domi- 
nated systems in the Universe. We show that, after successful rejection of interlopers, 
the line-of-sight velocity dispersions are substantially smaller than previous calcula- 
tions. With this, we use Jeans modelling to show that the DM content interior to 
300 pc appears to increase with satellite luminosity. Moreover, if we assume that 
MOND provides the true law of gravity, a correct identification of interlopers implies 
that these dSphs have reasonable values for their mass-to-light ratios, alleviating one 
of the claimed problems of MOND. 



t-H ; 1 INTRODUCTION 

, Impressive observations by I Walker et alj (|2007m and lMateo et al.1 (|200cf ) have accurately catalogued thousands of stars in the 
■ vicinity of the eight classical dwarf spheroidal galaxies (dSphs) of the Milky Way. By binning these line-of-sight (los) velocities 
as a function of projected radius and by computing the root mean squared (RMS) velocity within each bin, it is possible, 
through Jeans analysis, to reconstruct the mass profile with considerably less uncertainty than with only central los velocity 

dispersions. 

By fitting mass model s to these and to the ultra-faint dSphs (|Belokurov et al.ll2007l; iKoposov et al 
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120071 ; ISimon fc Gehal 



120071 ; IWillman et al.ll2006l ) that were discovered with the aid of the SDSS survev. lstrigari et alj (|2008f ) have proposed that 
' the mass within 300 pc of every dSph galaxy, over four orders of magnitude in luminosity, lies in a small range near 1O 7 M . 
k>( , This result is at odds wit h the naive expectat ion that more luminous galaxies reside in more massive halos; one therefore 
. needs to conclude that the lStrigari et al.1 l|200Sh result is of some fundamental significance to the dark matter particle which 
' creates the halos at high redshift, in which the dSph galaxies form, or to feedback processes which control star formation 
(|Koposov et al.|[2009h . 

Several other analyses have been made that rely on the los velocity dispersions, for instance the c alculation of t he dSph 
mass-to-light ratios M/L (the only truly fre e param eter) in Modified Newtonian dynamics (MOND) by [Angus (200^) and the 
mass modelling performed by I Walker" et all (|2009al ). 

An issue may exist with these studies because they depend crucially on the dSph (and all the stars within the projected 
tidal radius) being in virial equilibrium, except for some very obvious field stars which can and have been rejected. In general, 
there is no concrete way to know if a star is bound to the dSph or not. This is especially true of the cold dark matter 
(CDM) paradigm, since in principle we can have any density profile of dark matter halo which can simply be increase d to 
accommodate a fast star. Nevertheless, a re-analysis of the SDSS discovered dSphs Segue I l|Niederste-Ostholt et alj|2009l ) and 
Hercules (M. Wilkinson, private communication), that has been performed with alternative techniques to the one presented 
here, has demonstrated that neither of these two ultra- faint dSphs are as massive as previously thought, and unbound stars 
are frighteningly common. 

The point is that if the bound stars have RMS velocities of ~ lOkrns -1 and just a few unbound stars have velocities 
of (for instance) ~ 20kms _1 , the unbound stars dominate the combined RMS velocity. Since in Newtonian dynamics galaxy 
mass is proportional to the second power of the velocity dispersion (and to the fourth power in modified dynamics) it is 
absolutely critical to make sure the stars used in the Jeans analysis are bound. 

This is contrary to the mass modelling of spiral galaxies and clusters of galaxies. For instance, the masses of spiral galaxies 
are usually measured from the velocities of the neutral hydrogen gas which is forced to be on circular orbits and tidal features 
would leave very obvious signatures. For one thing, escaped gas along the los would have a totally different density to the gas 
inside the galaxy. Secondly, warps are relatively easy to identify in gas, but not from discontinuous stars. 
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The same is true for clusters of galaxies where the mass is inferred from the measured temperat ure and density of the 
X-ray emitting plas ma. Furthermore, cluster masses can b e measured with gravitational lensing (e.g.. [Clowe et al1l2006l ) or 
the caustic method ijDiaferio fc Gellerlll997l ; lDiaferidll999l l. which both ignore the state of equilibrium. 

Another important consideration is that the central regions of clusters are tidally influenced only slightly by their sur- 
roundings, both because of their mass dominance and the time scales involved. Similarly, most spirals exist blissfully in the 
field: consider even the "high density" of the local group spirals. At the opposite end of the spectrum, dSphs are in close 
proximity to their hosts (~ lOOkpc) and the tidal field is rapidly changing as their orbits are shorter (~ lOOMyr). Furthermore, 
the dSphs are extremely low surface brightness and extended, as opposed to the globular clusters, their counterparts in mass 
and position. 

As mentioned above, there is a method called the caustic technique which utilises the three known phase-space coordinates 
to delineate caustics in the los velocity-projected radius diagram which enclose the bound st ars. Although t he caustic technique 
was originally designed to study the mass profiles of galaxy clusters beyond the virial radius (Diaferio 2009), it can also function 
as an interloper identifier. This actually works independently of the gravity theory because the caustics are basically escape 
velocity curves and are generated in response to the proximity of stars to each other in the los velocity-projected radius 
diagram. Of course, in Newtonian gravity, this can be related to the enclosed mass, but it can also be used to confirm if stars 
are bound. 

This is of particular significance for the ultra-faint dSphs because lack of stars, greater susceptibility to tides and more 
severe impact of interlopers aggravate the problem, but in contrast to the ultra-faint dSphs, the classical dSphs have a highly 
statistically significant number of stars, which the caustic technique requires to operate effectively. 

After we have the confirmed stars as members, we can recompute the los velocity dispersion as a function of projected 
radius and perform our Jeans analysis in both MOND (to check whether the M/Ls are consistent with the stellar populations) 
and in Newtonian dynamics to investigate the claim that the mass within a radiu s of 300 pc is similar for every dSph. 

Here we make use of the recently published data collected by IWalker et all (|2009bh for Sculptor, Fornax, Carina and 
Sextans, as well as Leo I from lMateo et all l|2008l) . 



2 CAUSTIC METHOD 

The caustic technique (Diaferio 1999) was originally used to identify galaxy members and compute mass profiles in clusters 
of galaxies, from the central region to beyond the virial radius. This technique has proven to be very effective in estimating 
the mass profiles of galaxy clusters from N-body simulations, independently of the dynamical state of the cluster. Considering 
that the technique can identify members in potentially any system, we have modified it to be used for interloper removal in 
dSphs. The principle is the same, i.e., from the velocity diagram (los velocity vs. projected radius to the centre), the technique 
locates the curve which represents the escape velocity of the system. Therefore, the stars within these caustics are members 
of the system and they can be used to calculate the velocity dispersion profile. 

In the velocity diagram of an astrophysical system, the components populate a region whose amplitude decreases with 
increasing r. This amplitude A has been identified by Diaferio & Geller (1997) with the escape velocity along the line of sight, 
i.e. A 2 (r) = («e SC ios). To locate the caustics, it is necessary to determine a set of candidate member stars of the dSph. These 
stars will provide a centre, a size and a value for the los velocity dispersion, once they have been arranged in a binary tree 
(Figure [TJ, according with their pairwise energy 

3 Rp 2 mi + m J y ' w 

In this expression R p is the projected separation on the sky between the stars i and j, m^i are their masses, (Av) is their los 
velocity difference and G is the gravitational constant. 

In the following, we assign the same mass to every star. Therefore the tree automatically arranges the stars in potentially 
distinct groups with a single parameter only, the star mass. To get effectively distinct groups and to specifically define the set 
of candidates, we need to cut the tree at some level. This level sets the node from which the candidate members are hanging. 
The level is defined by the "cr plateau" extreme closer to the root, as described below. 



2.1 a plateau 

We can identify the main branch of the binary tree as the branch that emerges from the root and follows the nodes from 
which the largest number of leaves hang (highlighted in Figure [TJ . Following the main branch of the tree, from the root 
to the leaves, the los velocity dispersion crj os of the stars hanging from a given node x initially decreases rapidly. It then 
reaches a long plateau and again rapidly drops towards the end of the walk (Figure [5J ■ The plateau indicates the presence 
of a nearly isothermal system and the nodes x\ and X2 (where a\ os abruptly changes its slope) are good candidates for the 
system/substructure identification. The a plateau only appears when the stellar mass m in eq. ^ lies within a proper range 
[mi, mj], which depends on the dSph considered. In fact, each star is a tracer of the velocity and density fields, namely of the 
gravitational potential well of the system; therefore, the parameter m entering eq. {1} represents a fraction of the total mass 
of the system and is not necessarily close to the real mass of a star. The appearance of the a plateau guarantees that any 
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Figure 1. Binary tree of the stars in the field of Leo I. Stars are the leaves of the tree (at the bottom). The thick path highlights the 
main branch. The horizontal lines show the upper and lower thresholds used to cut the tree. Some nodes are labeled on the left side, 
with their number of descendants between parentheses. 
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Figure 2. Line-of-sight velocity dispersion of the "leaves" (stars) of each node along the main branch of the binary tree shown in Figure 
[T] Filled squares indicate the nodes x\ (left) and X2 (right). 



value of m within the range [7711,7712] is appropriate. By varying 771 within this range, we can estimate the systematic errors 
on the velocity dispersion derived from the caustic technique. Specifically, we apply the caustic technique for three values of 
m: m = mi, m = 7772, and a random value in between. 



2.2 Final members 



The candidate members are the stars hanging from the node x\. These stars determine the centre of the system: the systemic 
velocity is the median of the velocities, while the celestial coordinates of the centre are the coordinates of the two-dimensional 
density pea k. To find the p eak we compute the 2D density distribution on the sky f q (a,5) with the adaptive kernel method 
described in iDiaferiol (fl 999). It goes without saying that the caustic technique centres and the centres used bv lWalker et al.l 



and lMateo et all (|200St ) are indistinguishable. 



In the plane (r, v) (Figure O, where r is the projected distance and 7; the los velocity of each star from the dSph centre, 
the caustics are the curves satisfying the equation f q (r,v) = k. Here f q (r,v) is the star distribution in the plane, and k is the 
root of the equation 

(vL) K ,R = 4(v 2 ) . (2) 

The function (v 2 bc ) k ,r — f^ -^K.{ r ) l p{ r )dr / f Q ip(r)dr is the mean caustic amplitude within R, f(r) = J f q (r,v)dv, (v 2 ) 1 ^ 2 is 
the velocity dispersion of the candidate members and R is their mean projected separation from the centre. The stars within 
the upper and lower caustics are the dSph members, used to compute ai os (r). 
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Figure 3. Velocity diagram for Leo I (data from lMateo e t al. 2008). The black circles are all the stars in the catalogue, while the circles 
with a green star are the candidates selected by the binary tree. The red curves are the caustics identifying the members. The blue 
vertical line represents the mean separation R from the centre. 



3 RESULTS 

For the five dSphs, we separate the stars into radial bins of 100 pc (150 pc for Fornax). We then evaluate the standard 
deviation of the velocities of the star members in each bin around the systemic velocity, which gives us the los velocity 
dispersion (losVD) as a function of the projected radius. Figure U shows the velocity diagram of the five dSphs and Figure \S\ 
shows the cleane d losVDs for t he sta r members identified with the medium of the three sets of caustics we locate. The losVD 
as calculated by I Walker et~"aH \200ti ) are also overplotted for reference. 

Clearly, our losVDs are significantly smoother (in particular for Carina and Fornax) and lower (most noticeably Sextans 
and Sculptor) than the originals. The key point here is that the caustic technique has not significantly altered the magnitude 
of the losVD for the most distant and luminous dSphs (Leo I and Fornax). For these larger dSphs, it has merely pruned some 
outlying stars. Alternatively, for Sextans and Sculptor, the VDs have changed considerably and one can immediately see the 
non-symmetrical distribution of stars in the velocity diagram (Figure 2| before the caustic technique is applied. 

The nominal error bars of 0.5 km s -1 on each losVD bin are given as a guide only and as a reference for \ 2 statistics. In 
truth, the combined error from basic observational error on velocity, distance and using different star masses in the caustic 
technique (leading to different losVDs) is larger. Moreover, other errors orthogonally come from the dSph luminosity and 
light profile making a fully inclusive error analysis virtually impossible. What is certain is that if the caustic technique has 
successfully removed the majority of interlopers, the models reproducing the losVD across the projected radius should not 
systematically err by more than 0.5 km s -1 . 

To each of the five dSphs, we made a Jeans analysis to infer their DM halo properties or to inf er their M/L ratios in 
light of MOND. All the dSph parameters such as luminosity, distance and light profile are taken from I Angus! <|2009h . 



3.1 Dark halo properties 

To demonstrate how quantitatively different the DM halos required to fit the cleaned losVDs are, we us ed the cored dark 
halo density pro file p(r) =_£dm,o[1 + /tdm) 2 ] wm (cusped halos achieve just as good fits as shown in lAngus fc Diafericl 
2009) applied in lAngusl (120091 ) (where we refer interested readers to find a detailed description of the models) along with 
centrally isotropic, but radially decreasing, velocity anisotropics (/3(r) = f3 It - !")- The parameters of all the models are 

r +r /} 

given in Table 1. Fro m this we plot the enclosed mass within 300 pc for each of our 5 dSphs and compare with the analysis of 
IStrigari et all ([2008). One can see in Figure [6] that instead of a constant enclosed mass, independent of dSph luminosity, we 
find a clear trend for our sample of five. How significant this is will only be known once a similar analysis has been completed 
on the ultra-faint dSphs and the remaining 3 cla ssical dSphs. Sadly, some of the impact of this result is clouded because the 
original fits to the losVDs bv I Walker et al. (|2007l ) were made by marginalising t he cen tral data points within 300 pc, and this 



is why some of our data points still exceed the enclosed mass of IStrigari et"aU (|2008T ), even though our losVDs are lower. To 
counteract this, we overplot the enclosed mass from fits made by Angus (2009; where extra emphasis was placed on fitting 
the inner data points as well as the outer ones) and the trend for the enclosed mass to decrease after interloper removal is 
obvious. 

Clearly this might not be the end of the story, the other dSphs may remarkably align with 1O 7 M0. Furthermore, they 
may have different velocity anisotropies than presented here, although these are the simplest models that can achieve good 
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Figure 4. Velocity diagrams for our 5 dSphs with the 3 sets of caustics overplotted (see £|3.3| l. The dark pink diamonds indicate stars 
that are members of all three caustics, the black circles are stars that are members of no caustics and the blue stars are members of at 
least one, but not all the caustics. 

fits to the data. Interestingly, IWolf et all l|2009h used a technique which enables the measurement of the enclosed mass at a 
single radius (coincidentally the half-light radius) which is independent of velocity anisotropy. The half-light radius of Carina 
is close to 300 pc and their value is identical to ours. 



3.2 MOND 

In Milgrom's Modified Newtonian Dynamics (MOND) the only free parameter when Jeans modelling the dSphs is the mass- 
to- light ratio M/L of the stars. Naturally, this M/L must conform to population synthesis models with a reasonable initial 
mass function, which in the V-band is roughly 1-3 in solar units. Preliminary analysis of the dSphs by Angus (2008; whose 
method we precisely follow here) demonstrated that the more luminous dSphs had sensible M/Ls (around unity), but that 
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Figure 5. Linc-of-sight velocity dispersion profiles for the Walker et al. (2007) uncleaned samples (gray circles) and our cleaned samples 
(red squares). In the left-hand panels we overplot the best fit MOND models (parameters listed in Table 2) and in the right-hand panels 
we plot the best fit cored dark matter halo models (parameters in Table I). The two lines either side of the best fit are the best fit 
assuming the data points were increased (green) or decreased (blue) by 0.5kms — 1 . 
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Figure 6. The enclosed dark mass within 300 pc as a function of dSph luminosity. The purple dots are from Strigari et al. (2008) and 
the red squares are our cleaned samples presented in Figure [3] The empty circles correspond to the fits of Angus (2009) to the uncleaned 
Walker et al. (2007) data, using the same cored dark halos as we implement here. 
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Table 1. List of the parameters used in fitting the dark halos to the cleaned los velocity dispersions. All velocity anisotropics are 
parametrised as (B(r) = — ^ 0r ^ , except Sculptor which requires P(r) = f — r . In addition, we give the enclosed masses at 300 pc and 

r +r /3 r +r fj 

X 2 statistics for our fits. The first five rows refer to the caustics used in the main analysis, but the bottom four rows show the fits for 
two additional sets of caustics for both the Carina and Sextans dSphs, as discussed in £13.31 



the low luminosity dSphs, i.e. Carina, Sextans and Draco, had M/Ls that were too high. We have re-analysed the dSphs 
after removing the interloper stars and have found that it leaves all 5 dSphs with comfortable M/Ls (Figure [TJ. The velocity 
anisotropics used are parametrically given in Table 1. In particular, Sextans has dropped from a M/L of 9.2 to 2.5, and Carina 
from 5.6 to 2.0. Significantly, Leo I has not changed from 1.4 and Fornax has changed from 1.4 to 1.0. Therefore, it appears 
that unbound stars are a problem with a significance that increases the more nearby or low luminosity the dSphs are. 

An important point to grasp here is that the mere presence of the non-equilibrium stars is unusual in terms of the DM 
paradigm bec ause tidal forces are too fe eble at their current lo cations. There are numerical studies of tidal harassment of 
dSphs by e.g. iKlimentowski et all (|2008l ) and lLokas et"lil (2009), but all rely on the dSphs plungin g to pericentres that are 
a fac tor of five nearer the Milky Way (« 20 kpc) than expected from the observed p roper motions (jPiatek et al.l [20061 . 120071 . 
2008, Table Furthermore, given the observed d isk of satellites dKroupa et al. 2005) and that the Magellanic Clouds are not 
bound to the Milky Way (|Kallivavalil et aT]|2006bl lal) these orbits are improbable, so other mechanisms should be investigated. 

However, in MOND the dSphs are more loosely bound and the internal dy namics vary with distance from the Milky Way 
(due to the so-called external field effect, e.g. lMilgrom 1983, 1995; Angus 2008), thus giving a possible mechanism for plunging 
orbits. This idea is currently being studied with high-resolution N-body simulations in MOND (Angus et al. in prep). 

Furthermore, in the D M paradigm, it is not ob vious that stars can be easily stripped out of the dSph (although detailed 
studies are being made, e.g. IWetzel fc White! I2009T ). because the DM halo protects the dSph from tidal disruption, especially 
at the current distances of the dSphs. On the other hand, in MOND the extra gravity which binds the stars of the dSph 
together depends not only on the internal gravity of the dSph, but also on the external gravity exerted by the Milky Way. 
When this external gravity reaches a similar magnitude to the internal gravity of the dSph, the internal gravity decreases 
towards the simple Newtonian gravity of the stars. For this reason, dSphs are more susceptible to tidal disruption in MOND. 



8 A.L. Serra, G.W. Angus, A. Diaferio 



13 [ 



10 



Sculptor 



" "t 



10 



100 



Luminosity [10 5 L ] 



Figure 7. Inferred M/Ls in MOND using the old data (green asterisks, lAngudliOuih and our cleaned samples (blue squares). M/Ls 
consistent with population synthesis models exist between the red dashed lines. 
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Figure 8. A comparison between the los velocity dispersions found using different star masses in the caustic technique as discussed in 
93.31 The red data points corresponding to the best caustics are supplemented by blue and green points in each radius bin that signify 
the two other, extremal sets of caustics. A best fit model is fitted to each series of data points of the same colour. The top two figures 
are for the Carina dSph and the bottom two figures represent Sextans. The MOND fits are given in the left panels and the DM fits on 
the right. 



3.3 Star masses 

The star mass mentioned in 321 is the only free parameter included in the caustic technique. To estimate the response of the 
caustics to changes in the star mass, we produced three sets of caustics for each dSph by setting three different values of the 
star mass (as shown in Figure [8}. The highest and lowest values contain the range in which we find a recognisable a plateau. 
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Table 2. List of the parameters used in fitting the MOND models to the cleaned los velocity dispersions. All velocity anisotropics are 
parametrised by f)(r) = f3 tr~~2"> except Sculptor which requires /3(r) = j3 4 2 '^ r 5 . In addition, we give the M/Ls and \ 2 statistics 

for our fits. The first five rows refer to the caustic set used in the main analysis, but the bottom four rows show the fits for two additional 
caustic sets for both the Carina and Sextans dSphs, as discussed in i|3,3l 

This mass range is 1.3 - 1.4 x 10 7 M for Fornax; 2 x 10 6 - 2 x 10 7 M@ for Leo I; 1.8 x 10 5 - 1.4 x 10 6 M Q for Carina; 
3.5 x 10 5 - 7 x 10 5 A/ Q for Sculptor and 8 x 10 4 - 4.5 x 10 5 M o for Sextans. 

The differences in the three sets of caustics for each dSph are not significant (Figure and the number of stars lying 
within the widest caustics and outside the narrowest one do not cause a large variation in the losVD profiles. This emphasises 
that the method is robust. Nevertheless, to identify the spread in the context of the mass models, we fit models to the losVD 
from each caustic set for Carina and Sextans. The three models used for Carina and Sextans are barely distinguishable (see 
details in Table 2) and yield MOND M/Ls and DM enclosed masses that differ from the main caustics at roughly the same 
level as the errors we find by adding 0.5 km s -1 to the observed losVDs in the case of Sextans, but far lower for Carina. 



4 CONCLUSION 

Here we have app lied the caus ti c tec hnique to remove star interlopers in the 5 classical dSphs observed by IWalker et al.l 
( 2007. l2009al ) and iMateo et all (|2008[) . We have applied the technique to each dSph, taking into account the systematics 
introduced by the free parameter in the binary tree construction, namely the star mass. We have found that there are no 
major discrepancies between the caustics calculated from each star mass (as shown by Figure|?]and the two tables) in the mass 
range where a a plateau appears. Consequently, the velocity dispersion profiles are not strongly dependent on this parameter 
of the caustic technique. 

Using this technique we have shown that unbound/non-equilibrium stars are all too common in the low luminosity dSphs 
of Carina and Sextans. Removing these interlopers gives very different enclosed masses within 300 pc for fitted dark halo 
models. We have found that the mass increases by a fa ctor of three while the luminosity increases by a factor of ten. Thus, the 
common mass scale proposed bv lStrigari et al.1 (|2008h becomes unce rtain. Intrigui ngly, we also investigated the implications 
for MOND, which are altogether more positive. Original analysis by I Angus] l|200ct ) showed the M/Ls of the high luminosity 
dSphs were perfectly reasonable, but the low luminosity dSphs had abnormally high M/Ls, in particular Sextans and Draco, 
casting serious doubt on MOND's validity. After removing the interloper stars identified by the caustic t echnique, Sextans ha s 
a realistic M/L of 2.5 in solar units. We still await the final data release for Draco and Ursa Minor by I Walker et al.l l|2007l ). 
but eye inspection of their Figure 1 shows them to be likely highly contaminated by unbound stars. 
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